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We find that the evolution equation for the three-particle quark-gluon B-meson light-cone distri¬ 
bution amplitude (DA) of subleading twist is completely integrable in the large Nc limit and can be 
solved exactly. The lowest anomalous dimension is separated from the rest, continuous, spectrum 
by a finite gap. The corresponding eigenfunction coincides with the contribution of quark-gluon 
states to the two-particle DA so that the evolution equation for the latter is the same as for 

the leading-twist DA 4>+{uj) up to a constant shift in the anomalous dimension. Thus “genuine” 
three-particle states that belong to the continuous spectrum effectively decouple from (f- {uj) to the 
leading-order accuracy. In turn, the scale dependence of the full three-particle DA turns out to be 
nontrivial so that the contribution with the lowest anomalous dimension does not become leading 
at any scales. The results are illustrated on a simple model that can be used in studies of l/m^ 
corrections to heavy-meson decays in the framework of QCD factorization or light-cone sum rules. 
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I. INTRODUCTION 

B-meson light-cone distribution amplitudes (DAs) 
are the main nonperturbative input to the QCD de¬ 
scription of weak decays involving light hadrons in 
the final state Bi- In particular the leading-twist 
DA gives a dominant contribution in the heavy quark 
expansion and it received considerable attention al¬ 
ready il. Utility of the QCD factorization tech¬ 
niques depends, however, on the possibility to con¬ 
trol, or at least estimate, the corrections suppressed 
by powers of the &-quark mass that involve higher- 
twist DAs. This task is attracting increasing attention 
and in the last years there have been several efforts 
to combine light-cone sum rules with the expansion 
in terms of B-meson DAs [iMl- This technique al¬ 
lows one to tame infrared divergences which appear 
power-suppressed contributions in the purely pertur¬ 
bative framework and to calculate the so-called soft 
or end-point nonfactorizable contributions in terms 
of the DAs of increasing twist. One of the problems 
on this way is that higher-twist B-meson DAs involve 
contributions of multiparton states and are practically 
unknown. 

In this letter we point out that the structure of sub¬ 
leading twist DAs is simpler as compared to what one 
may assume from their general partonic decomposi¬ 
tion BUI- This structure is revealed by considering 
the scale dependence of the DAs in the limit of large 
number of colors, f oo, i.e. neglecting the 1/iV^ 
corrections to the renormalization group equations. It 
turns out that the evolution equation for the three- 
particle DA in this approximation is completely inte¬ 
grable and can be solved exactly. The lowest anoma¬ 
lous dimension is separated from the rest, continuous. 


spectrum by a finite gap. The corresponding eigen¬ 
function defines what can be called the “asymptotic” 
three-particle B-meson DA and has a relatively simple 
form. Most remarkably, it turns out that the higher- 
twist contribution to the two-particle B-meson DA 
4>- (uj) that is related to the three-particle DA by QCD 
equations of motion (EOM), is expressed entirely in 
terms of this “asymptotic” state, the states that be¬ 
long to the continuous spectrum do not contribute. 
As the result the DA </)_ (w) evolves autonomously and 
does not mix with “genuine” three-particle contribu¬ 
tions. The evolution equation for is the same 

as for the leading-twist DA (/>+(a;) up to a constant 
shift in the anomalous dimension. Finally, we discuss 
the evolution of the three-particle DA itself and its 
asymptotic behavior at small and large quark/gluon 
momenta which turns out to be nontrivial. This be¬ 
havior is illustrated on the example of a simple model 
that can be used in phenomenological applications. 

II. EVOLUTION EQUATIONS 

Following the established conventions Q we define 
the B-meson DAs as matrix elements of the renor¬ 
malized nonlocal operators built of an effective heavy 
quark field hy{0), a light (anti) quark and gluons at a 
light-like separation: 

= {0\q{nz)'^j5hy{0)\B{v)), 

= {Q\q{nz)i'y^hy{i))\B{v)) ( 1 ) 

and 

-2iF(^)$3(2;i,Z2,Ai) = 

= {Q\q{nzi)gGfj,,, {nz2)rF ( 0 ) 15 ( 1 ;)). ( 2 ) 
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Here is the heavy quark velocity, is the light¬ 
like vector, = 0 , such that n • w = 1 , F stands for 
an arbitrary Dirac structure, \B(v)) is the H-meson 
state, /r is the factorization scale and F{^) is the B- 
meson decay constant in the heavy quark effective 
theory (HQET). Wilson lines connecting the fields are 
not shown for brevity; they are always implied. 

The functions $+ and are the leading- and 
subleading-twist two-particle B-meson DAs [2, and 
$3 is the (lowest twist) three-particle DA that is the 
only one relevant for the present study. In notations 
of 14| $3 = A — 4'y. These three DAs are related 
by an EOM 


with respect to the SL(2) invariant scalar product [2lj 
($i|‘i> 2 ), = / V,z^Uz)<^2{z), (7) 

JC- 

where the integration goes over the complex coordi¬ 
nates z in the lower half-plane C_ : Im z < 0 and the 
integration measure is defined as 

Vjz = [i{z - z)Y^-^. 

TT 

Going over from quark/gluon coordinates to the cor¬ 
responding momenta 


(l + z5,)$_(z) 


$+(z)-|-2 f wdw^ 3 {z,w) (3) 

Jo 


ci>(z) = r 
^0 




( 8 ) 


that can be solved to obtain $_ as a sum of the so- 
called Wandzura-Wilczek (WW) term expressed in 
terms of $+ 0 , and a certain integral of the quark- 
gluon DA $ 3 . The latter contribution is nontrivial 
because it involves a function of two variables. We 
will demonstrate, however, that this complication is 
to a large extent illusory as the integral appearing in 
the EOM essentially decouples from “genuine” quark- 
gluon correlations. This simplification is exactly anal- 
ogous to what has been observed before [la - B^ for 
the structure function g 2 {x,Q^) in polarized deep- 
inelastic lepton-proton scattering. 

The following discussion is based on properties of 
the renormalization group equations for heavy-light 
operators under collinear conformal transformations. 
The corresponding generators read 


can be done easily making use of the following expres¬ 
sions 0 : 

= r(2j) S{k - k') , 

= r(2j) {ksf/^-^ J 2 ,-i{ 2 Vks). (9) 

Staying in coordinate space for the time being, we 
write the two-particle DAs as 

1 r°° ^ 

$+(z) = --/ 

•2^ JO 

poo ^ 

$_(z) = -- / (10) 

^ Jo 

and the three-particle DA 


S+ = + 2jz , = z^z +j, S- = -dz , (4) 

where j is the conformal spin, jq = 1 for the light 
quark and jg = 3/2 for the gluon. The genera¬ 
tors satisfy the standard SL(2) commutation relations 
[S'+,S'_] = 2S'o , [<5'o,5'±] = ±S± . We distinguish 

the generators acting on quark and gluon coordinates 
by the subscript Sq and Sg, respectively. 

The starting observation is that both the one-loop 
renormalization group equations (RGE) for the DAs 
and the EOM relations are invariant under special 
conformal transformations i,l3- It is therefore nat¬ 
ural to expand the DAs in terms of the eigenfunctions 
of the corresponding generator Q 

(z) = . (5) 


^3(21,22) = 


2_3 


ds I 


duuu^ 


zfz. 




1*2 


( 11 ) 


Here and below u = 1 — u. Inserting these expressions 
in the EOM relation ([3|) one derives for the expansion 
coefficients 


(/_(s,/r) = 0+(s, ^) — 2s^ / duuu(j)^{s^u, gJj. (12) 


Invariance under special conformal transformations 
implies that terms with different values of s cannot 
get mixed by the RGE. Thus the leading twist contri¬ 
butions 4)+{s,g) must have autonomous scale depen¬ 
dence: 


They form a complete orthonormal set 

{Q^P\Q^P)g=n‘^j)s^-^=5{s-s'), 

1 pOO - 

dss^^-^Q^p{z)Q^j\z') = 


( 6 ) 

(z — z'^j 


where [ 3 ,[^ 

£+{s, d) = ‘2Cf [In (/is) - V'(l) - 5 / 4 ] . (13) 
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The RGB for the three-particle DA $3 is more com¬ 
plicated, 

V ^ $ 3 ( 2 : 1 ,22, Ai) = 0, 

^ ^ (14) 

where the “Hamiltonian” Ti, to the one-loop accuracy 
is given by a sum of two-particle kernels 


and therefore the RGB (fHl) is “diagonal” in s. This 
symmetry alone is not sufficient, however, to find the 
solution since the problem has two degrees of freedom 
— the light-cone coordinates of the light quark and 
the gluon. It turns out, however, that for the leading 
contribution for a large number of colors 

n = NcM + n-Hm, (20) 


n = Hqg+ Hgh + Hgh ■ (l5) 

Bx plic it expressions for the kernels are known [^, [^ . 


[Hqhf]{zi) = - 


c KJO 


^ ■ 
a - 
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f{zi) - af{azi) 


[Hghf]iz2) = ^ 7 ( 22 ) - a^fiaz2) 




ln(*Ai22) - i /(22) !>, 
^ da r„ 


[Hqgy^Kzi,Z2) = f — 2ip{zi,Z2) - aipizf2,Z2) 
Uo « L 


-aV(2l,2^l)J - -ip{zi,Z2) 


~W< ' df3 P ip{z^2^ z^^), 


''C Jo 


(16) 


where 


^“2 = 0:2:1 -|- az2 , a = \ — a . ( 17 ) 

Note that in difference to [ 1 , [ 13 , I 13 | we include 
the QCD coupling in the definition of the quark- 
antiquark-gluon operator, G^i, :—>■ gGg,^. This redefi¬ 
nition affects the constant terms in the kernels. 

Bor our present purposes it is convenient to write 
these integral operators in terms of the generators of 
SL{2) transformations H, [13 


Hqg = N, 


(J,g + 3/2) + V' iJqg-^/2) - 2V^(1) - 3/4 


I 2 ( t\Joc- 3/2^ (Jqg 3/2) 

r(j„ + 3/2)’ 

Hgh = Nc [In {igS^) - V'(l) - 1 / 2 ] , 

1?,/. = -^ [In (*m5+) - ^(1) - 5/4] , (18) 


where Jqg is defined in terms of the corresponding 
quadratic Gasimir operator Jqg{Jqg — 1) = {Sq + SgY'. 
This representation makes manifest that the Hamil¬ 
tonian "H commutes with the generator of special con¬ 
formal transformations 


(19) 


there is an additional “hidden” symmetry. Namely, it 
is possible to construct one more “conserved charge”, 
Q 2 , that commutes both with Qi and the large- 
Hamiltonian H: 

[Qi,Q 2] = [Qi,H] = [Q2,H] =0. (21) 

Having two conserved charges for a problem with 
two degrees of freedom allows one to diagonalize the 
Hamiltonian, i.e. in our case find the multiplicatively 
renormalizable operators, without the need to solve 
the RGB equation (IT4l) explicitly. This property is 
known as complete integrability. 

The explicit expression for Q 2 can be found us¬ 
ing the formalism of the quantum inverse scattering 
method (QISM) [H ]: 


(22) 

In this approach the charges Qi, Q 2 appear in the ex¬ 
pansion of the element C(u) of the monodromy ma¬ 
trix for an open spin chain, C{u) oc Qi -I- Q 2 . The 
commutation relation [C(u) , H] = 0 can be verified by 
a direct calculation using the coordinate-space repre¬ 
sentation for the kernels as given in Bq. (HU) . or, more 
elegantly, with the help of the QISM techniques. This 
derivation will be given elsewhere . 

The “conserved charges” Qi, Q 2 and the “Hamil¬ 
tonian” H are self-adjoint operators with respect to 
the SL{2) scalar product ([7]): 

(V[/|(J)) = f ViZi [ 773 22 $*( 2 : 1 , Z2)$(2i. 2 : 2 ) . 

Jc- Jc- ^ 

It follows that they have real eigenvalues and can be 
diagonalized simultaneously: 

QlYs,xizi, Z2) = sTs.x( 2 i, 22 ) , 

< 02 ^ 8 , 2 ,( 21 , 22 ) = -Sx‘^Ys,xizi,Z2) , 
MYs^x{zi,Z2) = E{s,x)Ys^x:{zi,Z2) ■ (23) 

Note that we write the eigenvalues of Q 2 as a product 
—where s > 0 is an eigenvalue of Qi and x^ is 
a real number (but not necessarily positive), S R. 
This structure is motivated by QISM [l^l. The eigen¬ 
functions 1 ^ 8 , 2 ,( 21 , 22 ) are labeled by two “quantum 


Qi = *(5+ + ^;), 
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numbers”, s and x, and provide the basis of the so- 
called Sklyanin’s representation of Separated Vari¬ 
ables [ 2 ^. They can be found using the method de¬ 
veloped in [^ . 


ys,x(zi,Z2) = 


IS 


2_3 




zrz: 


1^2 Jo 


X 2 F 1 


2 2 


■ (24) 


The functions Yg^x are symmetric under reflection 
X —>■ —X. Since the eigenvalue x'^ has to be real, 
X can take real or imaginary values. It is possible 
to show that for imaginary x there exists only one 
normalizable solution corresponding to the particular 
value x = il2. For this special solution the hyper¬ 
geometric function disappears and the eigenfunction 
becomes very simple: 

„'„2 pi 

T®(zi, 02) =-r-T / . (25) 

zizi Jo 

This solution has the lowest energy 

Eq = E(s, cc = i/2) = ln(/rs)— '!/’(!) —1/4 (26) 


and can be interpreted as the ground state of the 
large-iVc Hamiltonian. It describes the “asymptotic” 
quark-gluon DA with the lowest anomalous dimen¬ 
sion. The state is normalized as 

(n^°^in®) = <J(«- 5 ')- ( 27 ) 

The eigenfunctions corresponding to real values of x 
belong to the continuous spectrum. They are orthog¬ 
onal to the ground state, = 0, and nor¬ 

malized as 


where i?+(s) is the anomalous dimension for the lead¬ 
ing twist DA $+, Eq. m, and A is a constant that 
does not depend on s: 

A = N^ + iV-i (ttVS - 3) . (32) 


The value of A coincides exactly with the gap be¬ 
tween the spectrum of anomalous dimensions of twist- 
three light quark-antiquark-gluon operators and the 
leading-twist quark-antiquark operators for a large 
number of derivatives, cf. [l8l| . 

A generic three-particle DA $ 3 (^ 1 , Z 2 , m) can be ex¬ 
panded in the eigenfunctions of the large-Ac Hamil¬ 
tonian 


^3{Zi,Z2,fl) = 


ds 


+ 


mis,d)Y}°Hz:i,Z2) 

^00 

dxr/(s,x,^)Yg^x(zi,Z2) (33) 


where the coefficient functions ? 7 o(s,/i) and r]{s,x,^) 
can be calculated using the SL{2) scalar product as 

Vo{s,fi) = , 

T]{s,x, fj.) = X ta.nhTTX (x'^+ ^j (Vs,a;|$3) • (34) 


They have autonomous scale dependence up to 1/Nj 
corrections: 


77o(s, fj.) = L^‘/^°R{s; /r, /tq) ?7o(s, do), 
r]{s, X, /r) = L'^^/^°R{s; /r, ^ 0 ) v{s, x, do) > (35) 


where L = as{fj,)/as(^io), /3o = - |n/. 


(V,.|nq.<) = 5(s - s')<5(3: - ■ (28) 

The corresponding eigenvalue (energy) is 

E(s, x) = \n{fis) + '0(3/2 -I- ix) + 0(3/2 — ix) 

-30(1)-5/4. (29) 

The gap between the ground state and the continuous 
spectrum 

AE = E(s, 0) - Eo = 20(3/2) - 0(2) - 0(1) (30) 

coincides with the gap in the spectrum of anoma¬ 
lous dimensions of twist-three light quark-antiquark- 
gluon operators with large number of derivatives, see 

Ref. Hi. 

The I/A/ corrections to the ground state energy 
£lo(s) = AcEo -I- l/Nc5Eo can be calculated in a stan¬ 
dard quantum-mechanical perturbation theory, evalu¬ 
ating the matrix element |(5]H1|V*^°)). The answer 
can be written as 

£o{s) = £+{s)+A + Oil/N!), (31) 


lx = Ac [ 0 ( 3/2 -I- ix) -I- 0 ( 3/2 - ix) + 27 b] (36) 


and 

R{s-fj,,no) = exp 


dr 


cuspi^^s ( t )) ln(Ts/so) 


'Mo 

2C p 


jY\ ^0 /mqs 
do) V So 


2Cp 

^0 


InL 




(37) 
is 


Here sq = e rc„,.s„(Q!.s) = ^Cp + 

the cusp anomalous dimension , and we used 


that 


F(p) = L-^/^°Fipio) . (38) 

Explicit solution of the evolution equation for the DA 
^ 3 (zi, Z 2 , d) in the large-Ac limit presents our main 
result. 

The two-particle DA $_ (z, pi) can now be recovered 
from the EOM relation (jS]). Using the expressions for 
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the eigenfunctions in Eqs. (HI, (HSl) one obtains for 
the relevant integral: 


/ wdw ^s{z,w, fi) = 
Jo 

1 

= - - sdse^’^/^ 



Vo{s,fJ.) 


1 f°° 

+ - / dxr]{s,x,fi) 2 F 1 
^ J — 00 



( 39 ) 



Remarkably, all terms involving the hypergeometric 
function vanish thanks to the identity 


du u 2 F 1 


— }y—ix,—\+ix 



( 40 ) 


that is related to the orthogonality condition 
= 0. Thus only the ground state (with 
the lowest anomalous dimension) contributes to the 
EOM relation ([3]). One finds 


r-z ^ noo 

J wdw^3{z,w,n) =J sdaF’^^^ T]o{s,ii), 


or, equivalently, 


duuu(j) 3 {s,u,fi) =-■^r]o{s,fi), ( 41 ) 


leading to the following very simple relation in the 
s-space: 


<(>-(s,m) = </'+(s,Ai)+?? o(s,m) . ( 42 ) 

Going over from quark/gluon coordinates to the cor¬ 
responding momenta 

poo 

^±(- 2 ^)= / duj e~^^^cj)±{uj), ( 43 ) 

PCO 

^3izi,Z2)= / 

Jo 

can be done easily making use of Eq. ([9]). The SL{2) 
scalar product in momentum space reads 


and the eigenfunctions of the evolution equation take 
the form 


b"s,x(wi,W2) — 


= — du^/Ui Ji (2^/uJiUs) UJ2 J 2 (2\/uJ2Us) 
Jo 


X\/sU2Fi 


-\-ix,-\ 


■ IX 


( 45 ) 


In this way one obtains the following expressions for 
the two-particle DAs in momentum space (cf. @, 0 ) 


4 >- (w,Ai) 


ds 4 >+{s, fi) \/uJs Jl{ 2 ^/cJs ), 

) ^ 

ds ^+{s,n) +r]o{s,^i) Jo( 2 \/ws). 

( 46 ) 


The scale-dependence of the coefficients (/)+(s,^) and 
770 (s,/i) differs by a simple overall factor 


</'-h(s, m) = R{s', Mo) <i}+{s, Mo), 
mis, m) = L^^^°R{s; M, Mo) Vo{s, Mo), ( 47 ) 


where i?(s;/x,^o) is defined in Eq. (1371) and A = 
Nc + 0{1/Nc) is a constant, see Eq. (I3T|) . In other 
words, the subleading twist contribution to 4)-(u},^) 
is suppressed at large scales as compared to the WW 
contribution by the universal factor that does 

not depend on the light quark momentum. To the 
0{1/N^) accuracy there is no mixing with “genuine” 
quark-gluon degrees of freedom. 

It is tempting to define the “asymptotic” quark- 
gluon DA $|®(2 i,02,m) S'® the contribution with the 
lowest anomalous dimension (for a given s): 


$^"(01,Z2,m) 



r]Q{s,^l)Y}^\zl,Z 2 ). ( 48 ) 


if, 9 ) =‘2 / —^/*(wi,W 2 ) 5 (wi,W 2 ) , ( 44 ) 

J Wi 


The corresponding expression in momentum space 
reads 


J 


())3’*(a;i,a;2,M) = ~‘^2\/o7r J ds^/sr]^^\s, fi) J (2-yiwiu) J 2 (2-v/sa^ 


UJ 1 UJ 2 
OJl + UJ2 


+ UJ 2 ) - /o(wi + W 2 ) 


-I- UJi 


/l(wi ■^U} 2 ) - /l(wi) , 


( 49 ) 
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where 


where 


foo 

fk{u})= dsT]^°\s,n){y/ujs) Jk{2y/uJs) . (50) 
^0 


S(^) 



s 


(53) 


III. ASYMPTOTIC BEHAVIOR AT SMALL 
AND LARGE MOMENTA 


The behavior of ^{z) at 2 —>■ 0 is determined by the 
small-s asymptotics of rjois). Assuming a power-law 

behavior 770 ( 5 ) s“ one obtains 


One of the reasons why the renormalization group 
evolution is interesting is that it gives insight in the 
expected behavior of the DAs at large and small mo¬ 
menta, which is important for the status of factoriza¬ 
tion theorems. Although one cannot make any rigor¬ 
ous statements on the shape of the DAs at low scales, 
it is usually assumed that the “true” DAs have the 
same asymptotic behavior as in perturbation theory. 
This assumption proved to be successful for modeling 
of parton distributions and DAs of light hadrons, so 
that it is natural to use the same logic for heavy-light 
systems. 

For small momenta there are no surprises. Using 
explicit expressions we find 

072,77) 

(071,072,77)0(a;|), (51) 

respectively. This behavior is in agreement with ar¬ 
guments based on quark-gluon duality M- If both 
quark and gluon momenta are small one obtains 

(?i^"(o7i,o72,7i )J dss'^r]^°\s,n). 

The large-momentum asymptotics is much more in¬ 
teresting. An inspection of the the asymptotic DA 
reveals that it does not decrease for large gluon 
momenta 072 —>■ 00 (because of the last term that is 
072-independent). As a consequence, integral over all 
momenta is ill-defined, and the normalization of the 
asymptotic DA to a matrix element of a local operator 
even at a single scale is not possible. This problem is 
seen even better in coordinate space. Using the defi¬ 
nition in (H51l and explicit expression for Ys^\zi, Z 2 ) 
one obtains 


•^Tizi,Z2) 


Z2 


Z 1 Z 2 

Z 1 -Z 2 



1{Z2) (52) 


^3{Zi,Z 2) ^ —dzidz^ZlZ2— - — ■ (54) 

22 2l - 22 

This function is not analytic at the origin {zi, Z 2 } = 0: 
the limit Zi —>■ 0 depends on the way the variables 
approach zero and exists only if a > 2. For a = 2 one 
gets 


$3(21,22) 2 [ 1 -|-Zi/z 2 ] -|-C>((zi/Z2)^,22) . ( 55 ) 

If the gluon coordinate 02 —t 0 and the quark position 
zi is kept constant, the singularity ^ I /22 cannot 
be avoided. It translates to the constant behavior at 
large gluon momentum, as seen explicitly from the 
momentum space representation. 

The singular behavior ^ I/ 22 , corresponding in 
physics terms to the instability due to gluon falling to 
the center of the color-Coulomb field, is not a special 
pathology of the asymptotic DA: the contributions 
of the continuum spectrum are even more singular, 
^ (I/^2)3/2±-, so that the corresponding momentum 
space DAs are increasing (and oscillating) functions 
of the gluon momentum. We are able to show that 
all such singularities are, however, spurious and can¬ 
cel in the sum of contributions of the asymptotic DA 
and the corrections. Most importantly, this cancel¬ 
lation is not spoiled by the evolution: The ~ I /22 
singularity is not generated at higher scales provided 
it is not present already in the nonperturbative ansatz 
at a reference low scale. This result implies that for 
small Z 2 , alias large 072 > /r, the hierarchy of con¬ 
tributions with increasing anomalous dimensions is 
lost; the leading large-a72 asymptotics of the “asymp¬ 
totic” DA is exactly cancelled by the contributions 
with larger anomalous dimensions. This pattern ap¬ 
pears to be unconventional, we are not aware of ex¬ 
amples of similar behavior for light quark systems. 

To this end we rewrite the expansion (1331) as the 
integral over imaginary axis x ^ —ix 


$: 


/•oo 1 riaa / q\ 

:( 2 i, 22 ) = J ds Y}°\zi,Z2){Yj°'^\^3) + — J dxx tanTTo; f --j Ys,2,(zi,Z2)(Ys,x|$3 ) 

1 l—l+ioc / g\ 

= TT / / dxxta,mTx ( - - Ys,x(2i,Z 2 )(Ys ^|$3), 

2% Jo J-i-ioo \ 4y 


( 56 ) 
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where 

'^s,xi^l,Z2) = Ys^_ix{zi,Z2) . (57) 

Note that the contribution of the asymptotic DA is 
taken into account in the second line of Eq. (15611 by 
moving the integration contour to the left of the singu¬ 
larity at X = 1/2, due to tauTra;. This representation 
remains valid after the scale dependence (1551) of the 


coefficients is taken into account: The anomalous di¬ 
mension 'j-ix (1361) is an analytic function in the strip 
—3/2 < Rea: < 3/2 and 7 (±i/ 2 ) gives the anomalous 
dimension of the asymptotic DA. 

In order to study the limit Z 2 —> 0 we write 

Ys,x(zi,Z2) = 'Ts^x{zi,Z 2) + Ts-x{zi, Z 2 ) (58) 

where 


Ts,a; (21,-22) = 


2.^3 


duuu^-^ 


r(2a 


z^z. 


1^2 -'0 


r(a: - i)r(ai + I) 




_1. _ 'T' ^ _ 'J' 

l-2x 


(59) 


and use that the coefficient (Ys^a:|$ 3 ) in (1561) is symmetric under x —x so that one can replace '¥s,x{zi, Z 2 ) —>■ 
2Ts^xizi, Z 2 ) without changing the value of the integral. For small Z 2 


Ts,x(2i,22) 


2i22 


Z 2 J 


r(2x)r(3/2 - x) 

r(-i -ha:)r(| -Ha;) L 


1 + 0 ( 22 ) 


(60) 


r 


so that the asymptotic behavior of the DA $ 3 ( 21 ,- 22 ) 
at -22 —>■ 0 is determined by the position of the closest 
to the origin singularity of the integrand on the real 
negative x axis, to the left of the integration contour 
at Rea: = —1. In this way the pole at x = —1/2 
corresponding the ^ 1/22 behavior is always avoided 
and the closest singularity appears to be at a; = —3/2, 
corresponding to $ 3 ( 21,22 —1^ 0) ^ const, unless the 
matrix element (Ys^ 2 ,|$ 3 ) is more singular. 


Taking into account the scale dependence amounts 
to the insertion of the RG factors (1551) under the in¬ 
tegral. In this way additional singularities appear 
corresponding to the poles of the anomalous dimen¬ 
sion (gni), 7 -ix = Nc[tp{3/2 -H x) -H ^(3/2 - x) -|- 27 ^]. 
The singularity closest to the origin is at x = —3/2 
so that if the initial condition for the evolution 
$ 3 ( 21 , 22 , H^o) corresponds to a constant behavior at 
22 —>■ 0 , it will be modified to $ 3 ( 21,22 —>■ 0 ) ^ 
ln(/j, 22 ), corresponding to a “tail” l/w 2 in momen¬ 
tum space. The same behavior was found previously 
for the leading twist DA It is easy to see that 

in the other limit $ 3(21 —>■ 0 , 22 ) ^ In(/X 2 i) as well, so 
that our final conclusion is that gluon emission gen¬ 
erates a radiative tail ~ l/wi and/or ^ I/w 2 of the 
three-particle DA ^ 3 ( 0 : 1 ,W 2 ,^) for both, large light 
quark and large gluon momenta. This is natural as 
the corresponding terms are present in the evolution 
kernels. The reason and consequences of such a be¬ 
havior have been discussed at length in the literature, 
e.g. so we do not need to repeat this discussion 
here. 


IV. A SIMPLE MODEL 


For the simplest phenomenologically acceptable 
model of the leading-twist B-meson DA at a low scale 
fi = Ho one usually takes Q 

(/+(w) = ^ ^+(s)=e-^^^ (61) 

where As is defined as 

^=r—</+(a;). (62) 

As Jo w 

The value of Xb is the most important nonper- 
turbative parameter in the QCD factorization ap¬ 
proach [3, i, with current estimates in the range 
Xb 'zz 300- 600 MeV @,[3. 

In the same spirit, we consider a simple model for 
the three-particle DA at a reference scale 

(j)3{u}i,u}2,Ho) = ^ wia;2 , 

^0 

1 / 3 ( 5 , u, ^ 0 ) = 7>3e"®‘^“ , (63) 


where ip 3 is a- constant that can be related to the 
matrix elements of local quark-gluon operators din 


7>3=g 


\2 


(64) 


The recent QCD sum rule calculation [T^ gives X% — 
X^jj = —0.03 ± 0.03 GeV^. The corresponding DA in 
coordinate space is 


_ 2^:3 _ 

(I -H iuJoZiYil -H iclIo22)3 


$3(21,22,/Xq) = 


(65) 














2 



Figure 1: The ratio /4^3 a function of uJ 2 /^o for 

several values of uji/ujq for the model in Eq. (|63|l . 
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the leading-twist DA. The higher-twist contribution 
^ y )3 in the s-space reads 

%(s,fio) =(67) 

and the DA at higher scales can easily be calcu¬ 
lated from Eq. 66]), gzl). The higher-twist contribu¬ 
tion is suppressed at large scales by an overall factor 
(as(/^)/Q!s(/ro))^/^° (I55|) as compared to the leading 
twist. 

Let us now have a closer look at the three-particle 
DA (IS51) itself. The asymptotic DA corresponding to 
this model is 




(wi — 2tLio) ^ 


03^1 

H-(w2 +uJi — 2ujo) 

LUq 


1 — e“o 


( 68 ) 


and the two-particle DA ^i)_(w, /xq) for this model can 
be obtained directly from the EOM: 


g—lii/As 

As 


2 

3a;3 


U} 

1 - 2 — 
Wo 


1^1 

2 Wq. 


,( 66 ) 


where the first term is the WW contribution related to 


The shape of (/)|®(wi, W 2 ) is qualitatively different from 
()) 3 (wi,W 2 ): it is not factorizable as a product of the 
distributions depending on wi and W 2 , does not de¬ 
crease at W 2 —t 00 and becomes negative for large wi. 
This is illustrated in Eig. 1 where we show the ra¬ 
tio (I>T/4>3 s-s a function of W 2 /W 0 for several different 
values of wi/wq. 

The quark-gluon DA at higher scales is given by 


00 

^ J dsR{s]fi,fio) j dxxtan7rx(^x‘^-^Ys^xiuJi,uj2)L^-'"^^°{YsM^3). (69) 

0 

where the a:—integration goes along the imaginary axis with Re(a;) = —1. For numerical evaluation it is 
convenient to consider first the corresponding function in the (s, u)-representation (HU, 


J dx X tamrx ^ 2 F 1 ^ ^ ^ 


and calculate the DA in momentum space as 




(70) 


poo nl ^ 

(/)(wi,W 2 ) = / dss / du (wiStt)^^^ (w 2 Su) Jl (2-y/MSWi) J2 (2y/usu!2) 4>3 {s, u,n)- (71) 

Jo Jo 


For our model ([63jl 

3 cos Tra: \ 4 / 

and the function (p 3 (s,u,^) has a factorized form 


where 


rirf ^ X / , TraJSmTTX / 2 f 2 1 

W(u,^i) = - / dx 3 - [X - - ] [X-- 

Jc cos2 TTX V 4y V 4 

-4 - X, -i -I- X 


U(j)3{s,u,fi) = (p 3 R{s-,fi,fio)e ( 73 ) 


X 2 F 1 


u 


J^l-ix/Po 

( 74 ) 
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U 


Figure 2: The function W(u,fi), Eq. (1741) . on three differ¬ 
ent scales /r = /ro = 1 GeV, /r = 2.5 GeV and ^ — 10 GeV. 

The function W (u, /i) oc u u, /r) at three different 
scales fj, = fiQ = 1 GeV, = 2.5 GeV and /r = 10 GeV 
is shown in Fig. [21 Note that u and u = 1 — u have 
the meaning of the momentum fractions carried by 
the quark and the gluon, respectively: 

UJl = UbJ , U}2 = UUJ , W = -h 0^2 , (75) 

so that the scale dependence visualized in Fig. |2] cor¬ 
responds to a redistribution of the total momentum of 
the light degrees of freedom such that at large scales 
the gluon carries a larger fraction. Note also that 
W(u,^) becomes slightly curved at m —)> 1 but still 
vanishes, or, equivalently, the DA diverges 

in the same limit, but the divergence is softer than a 
power 1/u. 

Finally, we show in Fig. [3] the (normalized) DA in 
momentum space 

^ ^2 

(j) 3 {u,UJ,n) =—(j) 3 {uuj,uuj,n) (76) 

^3 

as a function of the quark momentum fraction after 
the evolution to ^ = 2.5 GeV (solid curves) and at 
the initial scale fj, = fiQ = 1 GeV (dashed curves) for 
four different values of the total momentum ujIojq = 
{0.3,1,3,9}. 

One sees that for small total momentum the DA is 
rather strongly suppressed by the evolution whereas 
the shape is only weakly affected. On the contrary, 
for large momentum there is no suppression and the 


DA is strongly tilted towards small values of u, corre¬ 
sponding to small quark and large gluon momenta. It 
would be interesting to analyze the large-momentum 
behavior using the expansion of the type suggested 
in @ (see also 0). Such a study goes, however, be¬ 
yond the tasks of this paper. 

V. CONCLUSIONS 

To summarize, we have shown that the evolution 
equation for the three-particle quark-gluon B-meson 
light-cone DA of subleading twist is completely in- 
tegrable in the large Nc limit and can be solved ex¬ 
actly in analytic form. The most important result for 
phenomenology is that “genuine” three-particle con¬ 
tributions of quark-gluon states essentially decouple 
from the subleading-twist two-particle DA (^-(w) 0 
so that its properties are similar to the leading- 
twist DA. A similar simplification has been found 
before for the structure function g 2 {x,Q'^) in polar- 
ized deep-inelastic lepton-hadron scattering [l6l - ll^ . 
Based on this experience, we expect that “genuine” 
three-particle contributions do not contribute directly 
to many physical observables in B-decays at tree 
level because three-particle and two-particle twist- 
three contributions to the products of currents are 
typically related by Ward identities; hence they can¬ 
not have a different scale dependence. We also expect 
that a similar simplification of the renormalization- 
group dependence holds for twist-four distributions 
as well. This study is in progress [1^. In this way 
we hope to be able to identify important degrees of 
freedom in multiparticle quark-gluon distributions in 
heavy mesons that can be parametrized by a mini¬ 
mum number of nonperturbative parameters. This 
would present a step forward towards understand¬ 
ing of subleading corrections in powers of the heavy 
quark mass and ultimately allow one to increase sig- 
nihcantly the accuracy of QGD predictions for heavy 
meson (and baryon) decays based on the heavy-quark 
expansion and/or light-cone sum rules. 
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